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1 Introduction 

In this article we want to describe what is known about time independent 
spacetimes from a global point of view. The physical situations we want 
to treat are isolated bodies at rest or in uniform rotation in an otherwise 
empty universe. In such cases one expects the gravitational field to have 
no "independent degrees of freedom" . Very loosely speaking, the spacetime 
geometry should be uniquely determined by the matter content of the model 
under consideration. In a similar way, for a given matter model (such as 
that of a perfect fluid), there should be a one-to-one correspondence between 
Newtonian solutions and general relativistic ones. 

The plan of this paper is as follows. In Sect. 2.1 we collect the information 
afforded on one hand by the Killing equation obeyed by the vector field gen- 
erating the stationary isometry and, on the other hand, that by the Einstein 
field equations. Throughout this paper we assume this stationary isometry 
to be everywhere timelike. Thus ergoregions are excluded. We try as much 
as possible to write the resulting equations in terms of objects intrinsic to 
the space (henceforth simply called "the quotient space") obtained by quo- 
tienting spacetime by the action of the stationary isometry. Much of this is 
standard. But since none of the references known to us meets our specific 
purposes, we give a self-contained treatment starting from scratch. Since all 
the models we treat are axially symmetric, we add in Sect. 2.2 a second, axial 
Killing vector to the formalism of Sect. 2.1. 

In Sect. 2.3 we first introduce new dependent variables for the vacuum 
gravitational field, namely a conformally rescaled metric on the quotient space 
and two potentials, using which the field equations have an interpretation 
in terms of harmonic maps from the quotient space into the Poincare half 
plane. These potentials, originally due to Hansen, are used in our treatment 
of asymptotics in Sects. 3.1,2. We then formulate the boundary conditions at 
spatial infinity appropriate for isolated systems and prove two basic theorems 
due to Lichnerowicz on stationary solutions obeying these conditions. These 
theorems are manifestations of the above-mentioned principle concerning the 
lack of gravitational degrees of freedom. The first result, the "staticity theo- 
rem" , basically states that the gravitational field is static when the matter is 
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non-rotating. The second one, the "vacuum theorem" , states that spacetime 
is Minkowski when no matter is present. In Sect. 2.4 we restore c, the velocity 
of light, in the field equations and show that these tend to the Newtonian 
ones as c — > oo. 

In Chapter 3 we we study solutions only "near infinity". (Note that by 
the Lichnerowicz vacuum theorem, such solutions can not be extended to 
all of R 4 exept for flat spacetime.) In the Newtonian case such solutions 
are known to have a convergent expansion in negative powers of the radius 
where the coefficients are given by multipole moments. The relativistic situ- 
ation is slightly at variance with our statement at the beginning concerning 
the Newton-Einstein correspondence: namely, there are now, corresponding 
to the presence of two potentials rather than one, two infinite sequences of 
multipole moments, the "mass moments" which have a Newtonian analogue 
and the "angular momentum moments" which do not. One may now study 
the two potentials and the rescaled quotient space metric in increasing pow- 
ers of 1/r, where r is the radius corresponding to a specific coordinate gauge 
on the quotient space which has to be readjusted at each order in 1/r. 

The results one finds are sufficient for the existence of a chart in the one- 
point "compactification" of the quotient space (i.e. the union of the quotient 
space and the point-at-infinity), in terms of which yet another conformal 
rescaling of the 3-metric, together with a corresponding rescaling of the two 
potentials, admit regular extensions to the compactified space. As summa- 
rized in Sect. 3.2 one is then able to find field equations for these "unphysi- 
cal" variables which are regular at the point-at-infinity and in addition can be 
turned into an elliptic system. From this it follows that the unphysical quanti- 
ties are in fact analytic near infinity and this, in turn, implies convergence for 
a suitable 1/r-expansion (r being a "physical" radius) for the original phys- 
ical variables. Furthermore the structure of the unphysical equations yields 
the result that the (physical) spacetime metric is uniquely characterized by 
the two sets of multipole moments. 

It is remarkable that stationary vacuum solutions satisfying rather weak 
fall off conditions at spatial infinity, by the very nature of the field equations, 
have to have a convergent multipole expansion. We believe that the topic of 
far-field behaviour of time-independent gravitational fields is by now reason- 
ably well understood. The main open problem is to characterize an a priori 
given sequence of multipole moments for which the expansion converges. 

In Chapter 4 we review global rotating solutions. In Sect. 4.1 we outline 
a result due to Lindblom which shows that stationary rotating spacetimes 
with a one-component fluid source with phenomenological heat conduction 
and viscosity have to be axisymmetric. In Sect. 4.2 we describe a theorem of 
Hcilig which proves the existence of axisymmetric, rigidly rotating perfect- 
fluid spacetimes with polytropic equation of state, provided the parameters 
are sufficiently close to ones for a nonrotating Newtonian solution. In Sect. 4.3 
we present the solution of Neugebauer and Meinel representing a rigidly rotat- 
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ing infinitely thin disk of dust. In the final chapter we treat global nonrotating 
solutions. 

In Sect. 5.1 we outline the essentials of a relativistic theory of static elastic 
bodies. The remaining sections are devoted to spherical symmetry. It has 
long been conjectured that nonrotating perfect fluids are spherical whence 
Schwarzschild in their exterior region. In Sect. 5.2 we discuss the present 
status of this conjecture. A proof exists when the allowed equations of state 
are limited by a certain inequality. While this inequality covers many cases 
of physical interest, the Newtonian situation suggests that the conjecture 
is probably true without this restriction. In Sect. 5.3 we review spherically 
symmetric perfect fluid solutions. The final Sect. 5.4 gives a short description 
of self-gravitating Vlasov matter in the sperically symmetric case. 

In the subject of time-independent gravitational field of isolated bodies 
there are some topics we do not cover. We do not address the question of 
the conjectured non-existence of solutions with more than one body. (Miiller 
zum Hagcn J58| has some results on this in the static case.) Furthermore 
we limit ourselves to "standard matter" sources. Thus Black Holes are ex- 
cluded. (For this see the article of Maison in this volume.) We also could 
not cover the interesting case of soliton-like solutions for "non-linear matter 
sources" , starting with the discovery of the Bartnik-McKinnon solutions of 
the Einstein- Yang Mills system (see Bizon [pd|.) 
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2 Field Equations 
2.1 Generalities 

Let [M, g^v) be a 4-dimensional smooth connected manifold with Lorentz 

metric t^ t „ of signature ( V ++)■ We assume M to be chronological, i.e. 

to admit no closed timelike curves. Let £ M be an everywhere timelike Killing 
vector field with complete orbits. Thus we do not allow points where £ M turns 
null, i.e. we exclude horizons and ergospheres. It follows (see J26|) that the 
quotient of M by the isometry group generated by £ M is a Hausdorff manifold 
N and that M is a principal R 1 -bundlc over N. Furthermore this bundle is 
trivial, i.e. M is diffeomorphic to R 1 x N. The fact that this diffcomorphism 
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is non-natural (whereas of course the projection ir mapping M onto iV is) 
plays a role in the formalism we shall now develop. 

Let us introduce the differential geometric machinery necessary for writing 
the stationary Einstein equations in a way naturally adapted to As far 
as possible, we will be interested in quantities and equations intrinsic to N 
("dimensional reduction"). For a similar treatment see the Appendix of j23j. 
We define the fields V and uj\ v \ — Wu„x] by 

V := V < (2.1) 

^a:=3^V^A]- (2.2) 

The 3-form uj^x vanishes if and only if £ M is hypersurface orthogonal - in 
which case (M, g^) is called static. More important than uj^x will be the 
2-form t^ v , given by 

Given the fields a^ v and CcWa carry the same information, since 



crfj.u ■— ^^u\i X - (2.3) 



u^x = 3V- 1 ^a„x]- (2-4) 

Equ. (2.4) is obtained by expanding the identity Cl^i/Xp] = 0, which follows 
from (2.2), and contracting with In a similar way we obtain the relations 

^^x^ vX = ZV- x g^ v (2.5) 

co^x^ X =^lu pi/ xu puX ^- (2-6) 

We now invoke the Killing equation for £ p , i.e. 

A;flV = + V„f M - 0. (2.7) 

Expanding lo^x in terms of we easily see that 

V^ v = V- 1 [a flv + (V [tl V)^ ] ], (2.8) 

or, equivalently, 

<V = V 2 V [p (V- 1 ^ ] ). (2.9) 

In the static case we have — 0, whence there exist global cross sections 
given by t = const, where ^ = W^t. 
Equ. (2.9) implies that 

v^v- 2 ^) = o. (2.io) 

Clearly we have C^t = 0, where r is the 3-form given by t m „a = V^uj^x- 
By (2.10) and the identity £j = £Jdr + d(£Jr), this implies £\dr = 0. Since 
dr is a 4-form and £^ is nowhere zero, we infer in 4 dimensions that dr is 
zero, i.e. 

\7^(V- 2 u;„x P ]) =0. (2.11) 
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Equ.'s (2.10,11) are integrability conditions for the Killing equations (2.7) 
which are "purely geometric" in that they do not involve the Ricci (whence: 
energy-momentum) tensor. Now recall the relation 

V M V„& = -Ruxn'tp, (2.12) 
which follows from (2.7) and its corollary 

<?^V,V P ^ = - ( 2 - 13 ) 
From (2.2), (2.7) and (2.13) we find that 

VV = ^W R X]^, (2.14) 

which, using (2.8), implies 

V^V-V^) = 2V~%R X] ^ X - V- 3 *^*^, (2.15) 

where we have also used (2.6,7). Interpreting G^ u = Rp V — hg^vR as the 
energy-momentum tensor of matter, the r.h. side of Equ. (2.14) is zero iff the 
matter current, for an observer at rest relative to is zero. In that case, 
and provided that M is simply connected, there exists a scalar field u, called 
twist potential, such that 

^VA = ie MI ,A p V p a;, (2-16) 

and then (2.11) implies 

V' i {V~ 2 V ll u) = 0. (2.17) 

Note that, by virtue of £[ M w„Ap] = 0, w satisfies C^u = 0. 

Next, using the definition (2.1) and Equ. (2.8), it is straightforward to 
show that 

VpV^ = -2i? MA „ p £V + 2V~ 2 {a^a„ x - {V x V)^a v)x 

+ |W„W,V+^(VV) J ] (2.18) 

and whence 

V*V M V = -2i?^r + F-^W) 2 + 2V- 2 a„ u a^ u . (2.19) 

Now recall (see e.g. the Appendix of f23]| ) that there is a 1— 1 correspondence 
between tensor fields on M with vanishing Lie derivative with respect to £ M 
and such that all their contractions with £^ and ^ are zero - and ones 
of the same type on N. In the case of covariant tensor fields on TV, this 
correspondence is the same as pull-back under ir. Examples of such tensor 
fields on M are the scalar field V, the symmetric tensor field 

K v := g^ - V-%& (2.20) 
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and the 2-form a^ v = u) pvX £, X - Note that a pi/ can also be written as 

(r lu , = Vh/h v v 'V ll iZ v >. (2.21) 

The tensor is, of course, the natural Ricmannian metric on TV. The co- 
variant derivative £> M associated with h pL , acting, say, on a covector living 
on JV, is given by 

D li X v = h/h v "'V li ,X v ,. (2.22) 

Denoting by IZ^xa the curvature associated with Dp,, we find, using Equ. 
(2.9), that 

T^jxvXo = hfj 1 h v v h\ x h a a R^i v > y a > +2V~ 3 'a p V ox p -V~ 3 '(a X [pO v ] p -o p ^o „] A ) . 

(2.23) 

Since N is 3-dimensional, there holds 

°riy°\ P ] = 0, (2.24) 

so that 

n^ Xp = h/h/hx X 'hpp'R^ u ,x> P > + SV^a^axp. (2.25) 

Thus 

■Rp V = h/h/R Wv , - V-^R^Xp't"' t p ' + 3V- 3 a^x(T v x . (2.26) 
Using (2.18), Equ. (2.26) finally leads to 

- h/h/R^v, + hr-^D^v + 2v- 3 <j,x<j, x - \v- 2 (d p v)(d u v). 

(2.27) 

From (2.19) we deduce that 

D 2 V := hTDpDvV = -2R^? + ly-^DV) 2 + W^o^ . (2.28) 

We now make the following observation: when r M ... A is an arbitrary tensor on 
N, there holds 

h/ . . . fc A A 'Vv...v = (-VJ-^^Pl^V.i (2.29) 
Applying (2.29) to (2.15) it follows that 

(-Vy^D^K-V^a^} = hv'R^e- (2-30) 

Finally, projecting (2.10) down to N, it follows that 

D^(V- 2 a M )=0. (2.31) 

Given the spacetime (M, g pv ) with the Killing vector £ M , under the con- 
ditions stated at the beginning of this section, there are coordinates (t, x l ) 
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on M, such that the canonical projection tt takes the form tt : (t, x l ) i— > (x 1 ), 
with local coordinates on AT and such that the Killing vector £ M takes the 
form £ = d/dt. In terms of such coordinates tensor fields on N, say Ti...j(x), 
can be viewed as the tensor fields 

T li ... v {t,x)=S ll i ...6jT i ...j{x). (2.32) 

Since = V, there holds 

^dx^ = V(dt + tpidx*), (2.33) 

for some 1-form tpi. Note that, in the tangent space at each point (t, x l ) G M, 
the g MJ y-orthogonal complement of £ M is spanned by <pid/dt + d/dx l and the 
orthogonal complement of £ M in the cotangent space is spanned by dx % . From 
the definition h^ u — — V it follows that 

g^dx^dx" = V(dt + ip t dx 1 ) 2 + hijdx i dx j , (2.34) 

where V, <Pi, on the r.h. side of (2.34) are all independent of t. It is now 
straightforward to check that 

a^dx^dx" = 3{^V u C\])dx^dx l/ = V 2 d [l ip j] dx l dx j . (2.35) 

Thus cr M „, viewed as a tensor on N, is given by 

on = V 2 d[i<pj]. (2.36) 

In the static case t can be chosen so that ipi = 0. 

Conversely, let us start from the 3-manifold (N, hij , V, Oij ) with Ricman- 
nian metric h^, a negative scalar field V and the 2-form a^, subject to 

D {i (V- 2 a jk] ) = 0, (2.37) 

which corresponds to (2.31). Suppose, moreover, that N has trivial second 
cohomology. Then there exists a covector ifi on TV with 

ay =V 2 D [i i Pj] . (2.38) 

Define M = {t € R} x A~ and define on A~ the Lorentz metric g^ v by Equ. 
(2.34) and ^ by £ = d/dt. Then one checks that = V, that, under 

the projection tt : M — > N, is the pull-back of hij and that a^ u is the 
pull-back of (Tij = V 2 D[i(pj}. The fact that the product structure of M as 
M = R 1 x A^ is not natural is reflected in the above construction by the fact 
that <fi, solving (2.38), is given only up to <pi i— > ipi = ipi + DiF, with F a 
scalar field on N. Under this change <7 M „ given Equ. (2.34) remains unchanged 
only when we set 1 1— » t — t — F. 
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Given the fields (hij, V, ov, ) on N, we can define the fields r, r i: by the 
following equations: 

D 2 V = -2r+^V- 1 (DV) 2 + 2V- 2 <T l3 a i] (2.39) 

D'K-V)- 1 /^^ = (-V0 1/2 r, ( 2 -40) 
Hij = Uj + ^V^DiDjV + 2V- 3 a lk a 3 k - -^{D^D^). 

(2.41) 

It then follows from our previous considerations that the spacetime (M, 
satisfies 

R^dx^dx" = r(dt + f f dx e ) 2 + 2r i dx i (dt + (p e dx e ) + r ij dx i dx j . (2.42) 

In particular, iff r, ri, rij are all zero, (M, g^ u ) is a vacuum spacetime. In this 
case we refer to (2.39,40,41) as 'the vacuum equations'. 

For later use we record another form of the field equations 

G>„ - kT„„, (2.43) 

where 

Tpudx^dx" = r(dt + tpidx 1 ) 2 + 2n{dt + Lp j dx J )dx i + Tijdx i dx j , (2.44) 
and where we set 

g^dx^dx" = -e 2U (dt + Lp t dx 1 ) 2 + e^h^dx 1 dx\ (2.45) 

given by 

D 2 U =^(e- 4U T + f/)-e w u, ij u> ij (2.46) 

D l u) l3 = ne- iU T] (2.47) 
Hij = 2(D i U){D j U) - 2e 4U u) ik u) j k + h l3 e iu uo u u) kl + k(t 12 - h l3 f/). 

(2.48) 

Here 

u>ij = Uij = dfripj] (2.49) 
and indices are raised with . 



2.2 Axial Symmetry 

We now assume the existence of a second, spacelike Killing vector rf on 
(M,g^). There is the following identity 

4V^(r/ [p o VA] ) = -CnUuXp + GpR&ZxVp], (2-50) 
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where lo^ v \ is given by Equ. (2.2) and we have used (2.14). Suppose, in 
addition, that £ and r\ commute. Then the first term on the right in (2.50) 
vanishes so that 

4V(77 [p uv A] ) = 6eR^xV P ]- (2.51) 
In an analogous manner 

4V^ [P ^ A] ) = -tofR,*[vS\Vp], (2-52) 

where w'^ x is given in terms of r\ in the same way as w^x is given in terms 
of £. The r.h. sides of Equ.'s (2.51,52) are zero (at points where £ and rj are 
linearly independent) iff the timelike 2-plane spanned by £ and r\ is invariant 
under R^ u . These conditions will be satisfied when the energy momentum 
tensor is that of a rotating perfect fluid. We now assume that 77^ has an axis, 
i.e. vanishes on a timclikc 2-surface which is tangent to Then, and when 
the r.h. sides of (2.51,52) are zero, it follows that 

?7[ P <*VA] = i[p^ v \] = 0- (2.53) 

The relations (2.53), in turn, are nothing but the conditions for the 2-plane 
elements orthogonal to £ and 77 to be integrable ("surface transitivity of £ 
and 77"). The above result is due to Kundt and Trumper pS| ). 

For the purposes of Sect. 4.2 we need to transcribe the relations satisfied 
by 7y M on the quotient manifold N . Writing the 1-form rj^ = g^ v rj v as 

■q^dx^ = r)(dt + ip l dx 1 ) + t] i dx i , (2-54) 

so that 

^=^-^4t + ^ (2 - 55) 



the Killing equations 



are equivalent to 



V X dxg^ + 2 5A(A A)77 A = (2.56) 



^D.V = (2.57) 
2w y y = DiiV^ri) (2.58) 
Chij = 0, (2.59) 

where := Dutpj\. The surface transitivity conditions (2.53) get translated 
into 



f! {l D m = (2.60) 
riiiVjk] = 0. (2.61) 
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In particular, rf is a hypersurface-orthogonal Killing vector on (N, hij). Sup- 
pose, now, that the energy momentum tensor is that of a rigidly rotating 
perfect fluid, i.e. 

T^u = {p + P)u^u v + pg^ u , (2.62) 

with 

«m = + GVn), M = const, (2.63) 

and / is chosen so that is future-pointing and u^u^ = —1. With this 
specialization the quantities r, Tj, entering in the field equation (2.46,47,48) 
become 

t = f 2 ( P + p)(-e 2U + Q V f -pe 2U (2.64) 
n = f(-e 2U + n V )(p + p)nrj t (2.65) 
Tij = pe- 2U h l3 + f 2 Q 2 (p + p)r h r )j , (2.66) 

where r\i — h^rf . The normalization factor / is given by 

/ = [e- 2U {-e 2U + Q V ) 2 - r2 2 Wr 1/2 - (2.67) 

The field equations have to be supplemented by the Killing relations (2.57,58,59). 
Note that these imply that p and p are invariant under rf (in addition of 
course to being invariant under d/dt). Under these circumstances the con- 
tracted Bianchi identities, which imply that 

W„T^ V = 0, (2.68) 

boil down to the relation 

(p + p)f- 1 DJ = D lP , (2.69) 
the remaining condition, namely D J '{e~ AV tj) — 0, being identically satisfied. 



2.3 Asymptotic Flatness Lichnerowicz Theorems 

Before stating the conditions for stationary spacetimes to be asymptotically 
flat, we elaborate somewhat more on the vacuum field equations. First recall 
from (2.17) that, when M (or equivalently: N) is simply connected, there 
exists a field wonJV such that 

(Tij = \{-Vy/ 2 Eij k D k u, (2.70) 

where we have used (2.16) and 



e ijk dx i dx j dx k = {-V) 1,2 ^e^ Xa dx v dx x dx a . (2.71) 
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(Of course, the existence of oj could have also been inferred from (2.40) for 
ri — 0.) We now rewrite the vacuum equations in terms of the conformally 
rescaled metric hij (see (2.45)), given by 

hij = (-V)h ir (2.72) 

Then (2.39,40,41), together with (2.70) lead to 

D 2 V = V' 1 {DVf -V~ l {poS) 2 (2.73) 

D 2 lu = 2V- 1 (Dlu){DV) (2.74) 

and 

Kv = \v- 2 [(D l V)(D. J V) + (D^iD^l (2.75) 



or 

Gij = \v- 2 UpiVKDjV) + (D iU )(D jU ) - \ha\iDVf + (Doj) 2 \ 

(2.76) 

We can now give an interesting geometric interpretation of the vacuum equa- 
tions (2.73,74,76). Namely, let V be the Poincare half-plane with metric qAB 
given by 

q AB dz A dz B = V~ 2 (dV 2 + duo 2 ) (V > 0, -oo < u < oo). (2.77) 

Viewing (z 1 (x) , z 2 (x)) — (V(x),uj(x)) as a map from (N,hij) to (P, Qab)i 
one easily checks that Equ.'s (2.73,74) are exactly the conditions in order for 
this map to be harmonic, in other words 

D 2 z A + rg c z B tj z c \jh^ = 0, (2.78) 

where r^ c denotes the Christoffel symbols of qAB, composed with z c (x). 
The metric hij(x), of course, is not given, but has to satisfy (2.51). The r.h. 
side of Equ. (2.51), in turn, is nothing but the energy momentum tensor of 
the harmonic map. (V , qAB) can also be viewed as a spacelike hyperboloid in 
(2+l)-dimensional Minkowski space. Namely, define fields 

V 2 +L0 2 -1 UJ V 2 +UJ 2 + l 

$m = , $s = , $k = • (2.79) 

Then 

-<P 2 K +<P 2 I+ <P 2 =-l (2.80) 

Viewing {<Pk ,@m ,@s) as coordinates on R 3 with Lorentz metric 4(— d<fr 2 K + 
d$> 2 M + d$g), the induced metric under the map (2.79) is nothing but qAB- 
The fields #m,^s are the potentials first introduced by Hansen [^5| which 
we shall use in Sect. 3.1 and 3.2. 
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As with any harmonic map, we can associate a conserved current on 
N with any Killing vector on the target space V . Since V has 50(2,1) as 
isometry group, there are three independent such Killing vectors, namely 



I a d _ d 
dz A duj 



dz A dV duj 



We note in passing that the SO (2, 1) isometry of V is closely related to the 
"Ehlcrs transformation" discussed in the article by Maison in this volume. 
The conserved current ji associated with any Killing vector rj on V is given 

by 

ji = z A . lV B q AB . (2.84) 

Hence 

ji = V~ 2 DiLj (2.85) 

ji = V^Dy + r 2 wAw (2.86) 

ji = V^wDiV + (2V)- 2 (uj 2 - V^D.uj (2.87) 

i 

are all divergence-free on (N,hij). By (2.70) and (2.38), ji is also equal to 

ji=ei jk D m , (2.88) 
i 

and so the "charge" associated with ji is always zero. In the asymptotically 
flat we shall turn to later, (2.87) will be identically zero. 

The quantity (2.86) has the following spacetime interpretation (compare 
||23"1). Let S be a 2-surface in M which projects down to a smooth 2-surface 
on TV. Then there exist local coordinates (x^) = {t : x l ) such that S is given 

by 

*'V) = (0 > a;V)), A =1,2. (2.89) 

Now integrate the quantity e llupa V p £J y over S. After some computation one 
finds 

W(VT)|^ = {-V)-^{d iV - 2^y M ggJ (2.90) 
where, as before, — V^dutpjy Now, using (2.70), 

- 2(-VT 1/2 ^V + (-Vy^toD^ = D»'(we tf V). (2-91) 
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Thus, when £ is closed, the integral I of the expression (2.91) is given by 
(h = (-V)hij) 

I=[ {-V)- 1 {D l V + V- 1 ujD l cj)dS t . (2.92) 

The fact that this integral in vacuum only depends on the homology class of 
£ arises, in the spacetime picture, from the fact that 

V%f„] = 0, when R av = 0. (2.93) 

The quantity 

M = — I (2.94) 

is called the Komar mass of (M, guv)- For the Schwarzschild solution it coin- 
cides with the Schwarzschild mass when the "outward" orientation is chosen 
for dS\ 

We now come to the 



Boundary Conditions 

Recall that we require (M, g pv ) to be connected, simply connected and chrono- 
logical. Let, in addition, M contain a compact subset K and let M \ K be 
an "asymptotically flat end" . (The results of this subsection will remain to 
be true if M\K consists of finitely many asymptotic ends.) This means that 
M\K should be diffeomorphic to M R (R > 0) with 

Mr = {{x°,x l ) e R 1 x (R 3 \B(i?))} (2.95) 

with B(R) a closed ball of radius R. In terms of this diffeomorphism, the 
metric g pv in M \ K has to satisfy that there exists a constant C > such 
that (see @) 

M + \gn + r a \ 9flv - n MV \ + r 1+a \d a g^\ + r 2+a \d a d p g^\ < C (2.96) 

5oo < -C~\ g m < -C- 1 (2.97) 

VI'eR 3 g ij X i X j > C" 1 ^(X 1 ) 2 . (2.98) 

We assume a > 1/2. Furthermore we require R^ to be zero in M \ K. 
(This latter condition could be considerably relaxed.) It now follows that 
the level set x° = is a spacelike submanifold of M \ K which has a finite 
ADM- momentum (see ||). If p^ L is a timelike vector (which it will be 
for 'reasonable' matter except in the vacuum case), it now follows from the 
timelike character of that it has to be an asymptotic time translation, i.e. 



\C - A»\ + r\d a e i \ + r 2 \d p 8^\ < Cr- a , 



(2.99) 
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where the constants satisfy 

A^A^rjfiv < (2.100) 

(see §). Further more it follows from jjj, that in M \ K, or a subset thereof 
diffeomorphic to Mri for sufficiently large R' > R, there are coordinates 
(t, y l ) in terms of which g^ is again asymptotically flat with the same a > 1/2 
and so that £ M is of the form £ M d/dx^ — d/dt. Hence, in the coordinates 
(t,y l ), which we now call (t,x l ), the metric 

g^dx^dx" = V(dt + ^dx 1 ) 2 + h ij dx i dx j (2.101) 

satisfies 

IV+ll+rlSSjVI+^ISiajVI <Cr- a (2.102) 

M + r\dj<pi\ + r 2 \0 k d m \ < Cr~ a (2.103) 

\ha - <%| + rlflfcftyl + r 2 \d k d t h l} \ < Cr~ a (2.104) 

in M \ C. It follows that 

r\a l0 \+r 2 \d k a l3 \ < Cr~ a . (2.105) 

We remark that the time coordinate t, which is at first only defined in the 
open subset M\K of M, can be (Kobayashi-Nomizu |3(|) extended to a 
smooth global cross section of 7T : M — ► N. 

We now state and prove two uniqueness theorems due to Lichnerowicz 
which are basic for the theory of stationary solutions. 

Staticity theorem: Let (M, g Ml /,£ A ) be asymptotically flat with a > 1/2, £ M 
be an asymptotic time translation. If the matter is non-rotating relative to 
i.e. ri in Equ. (2.40) is zero, then the spacetime is static. 

Proof: From (2.40) we have 

D i [{-V)- 1 ' 2 a ij \=Q. (2.106) 

Contract Equ. (2.106) with ip> , using er^ = V^Dk^jI. It follows that 

!>'[(- Vj-VVfcV] = (-y)- 5 / 2 ^a ij . (2.107) 

Now integrate Equ. (2.107) over N. Since the term in brackets on the left 
is 0{r~ 2 ~ 2a ), the boundary term at infinity gives zero. Consequently cry = 
=> lo^x = 0. 

Remark: Since = 0, the field tpi is of the form ipi = DiF, where F = 
0(r 1_Q ). In the coordinates t = t — F, <? p „ takes the form 

g^dx^dx" = Vdt 2 + h l0 dx l dxK (2.108) 
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Vacuum theorem: Let (M, g^, £ A ) satisfy the conditions in the staticity theo- 
rem and let [M, g^) in addition be vacuum. Then (M, g M „) is the Minkowski 
space. 

Proof: Firstly, by the staticity theorem, we have that o~ij = 0. Using this in 
Equ. (2.39) for r = we have with v := (-V) 1/2 that 

D 2 v = 0. (2.109) 

By the maximum principle, or multiplying (2.109) by \x, integrating by parts 
and using /i — 1 = 0(r~ a ), difi = 0(r~ 1 ~ a ), we infer that /i = 1. Now 
Equ. (2.41) implies TZij = Hijkt = since dim N = 3. Since N is simply 
connected, it follows that (N, hij) is flat R 3 . Thus 

g^dx^dx" = -dt 2 + S ij dx i dx j . (2.110) 



2.4 Newtonian Limit 

Ehlers showed (unpublished, see j4i|) that one can write the field equation 
containing a parameter A = c~ 2 such that the equation remain meaningful 
for A = and then they are equivalent to the Newtonian equations. The 
variables for which this is true in the time dependent case have to be chosen 
in a quite sophisticated way. The stationary case can be treated in a direct 
and simple way as follows. 
We write the metric as 

g^dx^dx" = -er^ (cdt + ipidx 1 ) 2 + e^h ik dx i dx k (2.111) 

where we inserted "c" by dimensional analysis. The field equations decom- 
posed in section 2.1 

R M u = ^§-{T*» \ T 9»») (2.H2) 
for the energy momentum tensor 

r M „ = c 2 r(cdt + ifidx 1 ) 2 + 2cTi{cdt + Lp j dx j )dx i + T io dx l dx j (2.113) 
become 

D 2 U = 4mG{e~^T + c _2 t/) - u>ijU) ij (2.114) 

Z/wy- = S-rrGc^e'^Tj (2.115) 

Rij = 2c~ i D i UD J U - 2e^~ ui ik uij k + li i3 e^ ui k iCj kl + 87rGc~ 4 (ry - hijf{) 

(2.116) 

Considered as equations for U, hij, <pi, r, 7$, , the equations (2.114-116) have 
a limit for c — > oo. 

In the static case the limit is 



D 2 U = AttGt (2.117) 
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R lk =0 (2.118) 

Hence we obtain immediately that the metric hik of the quotient is flat and 
therefore (2.117) is the Poisson equation of Newton's theory. The connection 
has also a limit and the only non vanishing Christoffcl symbol is r% t — D l U . 
The equation of motion V = becomes in the limit the Newtonian 
equilibrium condition I) 3 Tij = —tDjU. 

Now to the stationary case: Because the right hand side of (2.115) vanishes 
in the Newtonian limit, the Lichnerowicz theorem implies that u>ik — which 
in turn implies by (2.116) that Rik = whence the metric on the quotient is 
again flat. 

For the metric written in the form (1) the connection has no limit for 
c — > oo. If we use however as a consequence of the field equations that Ujj = 0, 
the connection has a limit and the equations of motion become the Newtonian 
equilibrium conditions, i.e. 

DWi = , D j Tij = -tD 3 U . (2.119) 
2.5 Existence Issues and the Newtonian Limit 

The fact that the equations can be written to contain A = c~ 2 in such a 
way that they are analytic in A and are the Newtonian equations for A = 0, 
suggests to use this structure for existence theory. In this section we will make 
some remarks about the static case. In section 4.2 an existence theorem for a 
rigidly rotating body by Heilig will be discussed which exploits the fact that 
the equations have a nice Newtonian limit. 

To obtain partial differential equations for which there is an existence 
theory we write (2.104) and (2.105) in the static case in harmonic coordinates 
on N, defined by D 2 x l = 0, for the unknowns U and Z %d defined by h lJ — 
S ij + \ 2 Z^ and obtain: 

AU := 5 l] d l d 1 U = AttGt + A(X, r, Tij , Z ij ) (2.120) 
AZ 11 = -Ad k UdiU5 tk 5 l3 - 167rG(f 4J - fS ij ) 

+ \ 2 B^(X, r, r kl , Z k \ d m Z kl , d m d n Z kl ) (2.121) 

here we used the well known expression for the Ricci tensor in harmonic 
coordinates 

fi*3 = _\^l dkdi y i ij + H ij( d ^ Qfy (2.122) 

where is quadratic in the first derivatives of h lJ . As usual we call (2.120), 
(2.121) the reduced field equations. These form a quasilinear elliptic system 
with the property that for given small sources r, Tij of compact support and 
small A there exist unique solutions U, Z % ^ which tend to at infinity. 

In particular we can choose for r, Tij a Newtonian solution and determine 
then for small A a relativistic solution of the reduced field equations which 
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have a Newtonian limit. Is is to be expected that the solution will be analytic 
in A. Then the Taylor expansion in A can be considered as a converging post 
Newtonian expansion. 

A solution of the reduced field equations is only solution of the field equa- 
tion if it satisfies the harmonicity condition or cquivalently if V^T^ = 
holds. To solve the reduced equations and the equation of motion is a much 
harder problem. It makes only sense once a matter model is chosen. 

In the static case only matter models of elasticity lead to new interesting 
problems because, as we will see in sections 5.2 , static fluids are spherically 
symmetric and can be investigated by ordinary differential equations (see 
Sect. 5.3). Some remarks on static, small self gravitating bodies can be found 
in Sect. 5.1. 

For a stationary rigidly rotating fluid Heilig has given an existence the- 
orem by perturbing away from a Newtonian solution. We will describe this 
result in Sect. 4.2. 

3 Far Fields 

3.1 Far-Field Expansions 

While, as we have seen, little is known so far about globally regular, asymp- 
totically flat solutions to the stationary field equations with reasonable matter 
sources, there is an almost complete understanding of the behaviour of gen- 
eral asymptotically flat solutions near spatial infinity, which we now describe. 
Here the quotient manifold N is of the form 



On N there are given the fields (hij, V, ipi) satisfying (2.39,40,41). The whole 
discussion is "local-at-infinity" . In particular one has to allow for R in B(R) 
to be made suitably large, as one proceeds. We will do so tacitly without 
changing the letter "i?" . The Einstein equations are given by 



By the asymptotic conditions (2.102-104), u> tends to a constant at infinity. 
Subtracting this from lo, and calling the result again lo, we find that 



N = R 3 \ B(R). 



D 2 V 
D 2 cu 





uj\ + r\diU)\ + r 2 \didjtv\ < Cr 



— a 



(3.4) 



In short, we have that 



V = -l + 0(r- a ), u = 0(r- a ), hij-dij =0{r- a ), 1 > a > 1/2 

(3.5) 
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and that these relations may be differentiated twice. The condition a > 1/2 
could be relaxed (see Kennefick and O Murchadha Qj). It now follows that 

AV = 0{r- 2 - 2a ) (3.6) 
Auj = 0(r- 2 - 2a ). (3.7) 

Since the r.h. sides of (3.6,7) decay stronger than 0(r~ 3 ), it follows from 
standard results in potential theory p3|j73] | that there exist constants M, S 
such that 

2M 

V = -l + +0(r- 1 - a ) (3.8) 

r 

o,c 

lo = — + 0(r- 1 -°). (3.9) 
r 

But the existence of tpi in Equ. (2.36) implies that S has to be zero. 

The equation (3.3), which involves second derivatives of the metric h^, 
yields 

A(ka - ^5ijk) - 2d (i r 3) + Sijde^ = 0(r- 2 - 2a ), (3.10) 
where fcjj = hij — 5ij , k = ha and 

r t = d,h l3 - ^dik. (3.11) 

This equation can be viewed in two ways, both of which recur in the higher- 
order steps leading to the theorem below. Firstly, in the gauge where Fi = 0, 
i.e. the linearized harmonic gauge for hij, it is an elliptic equation, namely 
essentially the componentwise Laplace equation, for the leading-order part of 
kij. Secondly, Equ. (3.10) can be rewritten as 

£i£ m £jn P did n k mp = 0(r~ 2 ~ 2a ), (3.12) 

which expresses the fact that the linearized Riemann tensor of hij decays 
faster than 0(r~ 3 ). Note that (3.12) makes essential use of the three-dimensionality 
of space. It follows [f73j that there exists gt = (9(r 1_Q ) such that 

= d l9j + j9l + 0{r- 2a ). (3.13) 

Thus the leading-order contribution to the metric hij is "pure gauge" . 
To next order in 1/r one finds that there is a gauge, namely 

r t = 0(r- 3 - a ), (3.14) 

for which 

V = 1 + ~ $~ + — + Oir- 1 - 2 *) (3.15) 

w =^- + 0(r- 1 - 2a ) (3.16) 
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and for which hij can be brought into the form 

h 3 = _ M2 ^-^) + ofr- 1 -**). (3.17) 

In the above M, Mi, Si are constants. All indices are lowered and raised with 
Sij. When M ^ one can, by a rigid translation, arrange for M, = 0. In that 
case the metric <? M „ obtained from (3.15-17) coincides, to order 1/r 2 , with 
that of the Kerr spacetime with \S\ = —Ma, M being the mass and a being 
the Kerr parameter. 

In order to extend the above result to higher orders in 1/r, it is convenient 
to replace (V, ui) by some other choice of scalar potentials. One choice, due 
to Hansen ^5j, is to set (see Equ. (2.79) 

V 2 +cj 2 -1 , s 

4>m = ^ (3.18) 

<t>s = -£ (3.19) 



V 2 +uj 2 + 1 

i>K = - 



(3.20) 



W 

It then turns out that (<p a ) = {4>m ,4>Si4>k) ah satisfy 

D 2 <j) a = 2K<j) a . (3.21) 
Then one has j73| the following 

Theorem: There exists a gauge, namely that where = 0(r~ m ~ 1 ~ a ), for 
which there are constants A ... ,B ....... ,G .. . such that 

i XJ %\—%l dj ••--Ij n oo/ _ m +i_2a\ fo oo\ 

<t> M = <L Tumi + (r ) (3.22) 

e=o 

ZT 1 ^1 f'H 

^ - <L + O (r ) (3.23) 

= | + E " / V" r " + 0°°(r-+ 1 - 2 «) (3.24) 

£=0 



Note that £ = M 

m 

hij = + 



£=2 



+ 0°°(r _m+1_2Q )J 
(3.25) 
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All constants are symmetric in their i\ . . . indices. The constants D are also 
symmetric in i and j. The constants C . . . appear only for m > 3, the con- 
stants D . . . only for m > 4. The symbol 0°°(r k ) means that the quantity 
in question is of 0(r k ), its derivative is 0(r fc_1 ), a.s.o. Furthermore all con- 
stants are determined by the tracefree parts of E . . . , F . . . in a way which 
does not depend on the solution at hand. The tracefree parts of E . . . are the 
analogues of the Newtonian multipole moments. The constants F . . . play 
an analogous role for the "angular-momentum aspect" , which does not have 
a Newtonian counterpart. The three- metric hij, for reasons explained after 
(3.14), has no independent degrees of freedom. 

This theorem shows, in essence, that any stationary, asymptotically flat 
solution to the Einstein vacuum equations is uniquely determined by the 
"moments" E...,F.... However no statement concerning convergence of 
series like the ones appearing in (3.22-25) can be made. In order to do that 
it is necessary to use "conformal compactification" of three-space N. 



3.2 Conformal Treatment of Infinity, Multipole Moments 

Before turning to the situation G.R., it is instructive to recall the Newtonian 
situation. Suppose we are given a Newtonian potential near infinity, i.e. a 
function <f> with 

A<j) = onR 3 \B(R). (3.26) 
Extending <j> smoothly to all of R 3 , we thus have that 

A(j> = <iirp withpeC oo (R 3 ) (3.27) 

and cj> — > at infinity. Thus (j> is of the form 

<t>{*) = - I r^\ dx '- ( 3 - 28 ) 

J R 3 \X-X'\ 

In R 3 \ B(R) this can (see e.g. pj|) be expanded in a standard fashion in 
powers of 1/r. One obtains an expansion of the form 

with E i± i £ totally symmetric and tracefree. One shows ]35| , |73| that this series 
converges absolutely and uniformly in R 3 \ B(R) for sufficiently large R. 

As a warm-up for G.R. it is useful to rephrase the Newtonian situation 
using "conformal compactification" . First observe that there is a positive 
smooth function fl on N — R 3 \ B(R) with the following properties. The 
metric 

hij = f2 2 S lj (3.30) 
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extends to a smooth metric on the one-point compactification 

N = N(J {r = 00} = NU {A}, (3.31) 

where 

f2\ A = 0, Ai?U = (3.32) 

and 

DiDjD - 2hij = 0. (3.33) 
To prove this, take Q — 1/r 2 and introduce 



x 



(3.34) 



as coordinates on N. One also sees that hij is again the standard flat metric 
in the coordinates x % . (This would also follow from (3.33) and the standard 
formula for the behaviour of Rij under conformal rescalings.) As for the 
potential, rewrite (3.26) as 



l2 n 



Z^--U = Q, (3.35) 
and observe that the operator in (3.35) obeys 

[d 2 - =Q-^(d*- (3.36) 

when hij = Q 2 hij and <j) = fi^ 1 ^ 2 ^ for arbitrary Q > 0. Thus we again have 

(5 2 -^ = 5 2 = O, (3.37) 

at first only on N. 

In the case of G.R. we were unable to prove convergence of the multipolc 
series, but only an asymptotic estimate like 

= E Eil "f2;i XU + 0°°(r-^- 2 "). (3.38) 



i=0 



But, from (3.38) for m = 4, it follows immediately that (j> extends to a C 3 - 
function on TV. Thus, by continuity 

D 2 4> = Q on N. (3.39) 

But it is a standard fact that solutions to the Laplace equation and, more 
generally, for nonlinear elliptic systems with analytic coefficients [d5[, are 
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themselves analytic. Thus 4> has a convergent Taylor expansion at the point 
A. But this is nothing but (3.29) in inverted coordinates. Furthermore the 
multipole moments Ei x .„i t can now be viewed as the Taylor coefficients of <p 
at A. It follows from (3.39) that they have to be tracefree, and it is trivial 
that they determine <p uniquely. 

Suppose fi is just required to satisfy (3.32,33). Then, given hij, there is 
in (hij , fi) the following 3-parameter gauge freedom 



fi' = u>fi, 
hit a = w 2 h 



(3.40) 
(3.41) 



where 



u) = (l-WDifi + b l b l fi)- 



(3.42) 



with Dib 3 = 0, which, in the compactified picture, corresponds to the freedom 
of choosing an origin in the "physical" space R 3 , w.r. to which the inversion 
x 1 = x l /r 2 can be made. Therefore the Taylor coefficients of U at A behave 
under (3.40,41) in a way which precisely corresponds to their dependence on 
the choice of origin. 

In G.R. it is impossible to require a conformal compactification for which 
(3.33) holds everywhere. We call a 3-metric hij on a manifold N = R 3 \£?(i?) 
conformally C , when there exists a C fc -function fi > on N such that 



fi 2 hij extends to a C k -metric on N 



N U {A} and 



fi\ 



DM 



(3.43) 



(DiDjfi - 2hij) 



= 0. 



(3.44) 



A scalar potential (f> is called conformally C k , when 4> = fi^ 1 / 2 ^ extends to a 
C fe -function on N. Given (3.43,44) there is now a much larger gauge freedom 
involved in constructing the unphysical from the physical quantities, namely 



fi' = wfi, 



hL 



LO hi 



-1/2, 



(3.45) 



where u> satisfies ui\a = 1- Now consider, following Geroch p2| ], this recur- 
sively defined set of tensor fields on N: 



Pa 
Pi 



!l...l m +l 



= Di4 
= TS 

TS 



DiDi 



D, 



P 

l+1 r i 1 ...i„ 



s(2s-l) 



V P 



(3.46) 
(3.47) 

(3.48) 
(3.49) 
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where TS denotes the operation of taking the symmetric, trace-free part. It 
turns out that the tensors 



(3.50) 



behave under (3.45) in exactly the same way as the Newtonian moments 
under the restricted gauge freedom (3.40-42) with 6j = Dhjj\a- Thus the Ricci 
terms in (3.46-49) cancel out unwanted dependencies from higher-than-first 
derivatives of ui at A. 

Now return to the expansions (3.22-25) for some fixed m > 1. Performing, 
again, an inversion x % = x l /r 2 and setting, in these coordinates, 



'M ■ 



4>s = ft 1/2 4>s, 



(3.51) 
(3.52) 



with fl = 1/r 2 we find that (0m, 4>Si hij) are all C m . Furthermore ft is C°°. 
Thus we have obtained a C m conformal compactification. Our proof would 
be complete if we could find an elliptic system satisfied by (hij, 4>m, <fis) or 
quantities derived from them. Doing this is not completely trivial. We explain 
the essentials in the static case where 4>s = 0. Thus 



V 



D 4> M = 2TZ(j>M 
2 



1 + 40 



^-(Di(j) M )(Dj<j)M)- 

M 



(3.53) 
(3.54) 



Let us assume that M ^ 0. Define, instead of 1/r 2 as above, a conformal 
factor also called fl by 

[(-V) 1 ' 2 - l] 2 



(3.55) 



(-V)V 2 ' 

It is not hard to sec from (3.22-25) that this yields a C m -compactification 
(4>M,hij) where, however, we have for convenience replaced (3.44) by 



2 ~ 

D i D i Q -^2 h H 



o. 



(3.56) 



It is useful to employ, as the scalar variable in the unphysical picture neither 
4>m nor f2, but the quantity a defined by 



a := 



(-V) 



1/2 



, 2 



_(_y)i/2 + 1 
After some labor we find from (3.53,54) that 



(3.57) 



K = 



(3.58) 
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and 



- <r(l - a)Tl v} = D t D 3 a - ]-h i:j D 2 a. (3.59) 



The scalar a satisfies 

tr\ A = 0, D it r\ A = 0, V 2 <t\a = ^. (3.60) 
Taking a "curl" of Equ. (3.59) we obtain 

(1 - cj)D {l K ]]k = 2{D {i a)R j]k - h k[l R j]e D e a. (3.61) 

If we take D l of the quantity D^{R,^ k and use the Ricci and Bianchi identities 
we find the relation 

(3.62) 

Using that 1Z is zero and Equ. (3.62), writing IZij = r^ , and using (3.59) to 
eliminate second derivatives of a, we obtain an equation of the form 

D 2 Tij — nonlinear terms, (3.63) 

where these nonlinear terms depend at most on , a and their first deriva- 
tives and on D 2 a. We call D 2 a = p. From the divergence of Equ. (3.59) we 
infer that 

pa = \{Do) 2 , (3.64) 

and from this after some work that 

D 2 p = 3ct(1 - i r R./R^ + 37^(5V)(Z)V). (3.65) 

Now Equ. (3.63) can be completed as follows: 

TZij = Tij (3.66) 

D 2 Tij = nonlinear terms (3.67) 

D 2 a = p (3.68) 

D 2 p = 3<r(l - <7 2 K, ! K , > + miji&a^&a). (3.69) 

Going over to harmonic coordinates, the "non-elliptic" terms in the expres- 
sion of IZij in (3.66) in terms of the metric go away, and the whole set of 
Equ.'s (3.66-69) becomes an elliptic system. Note that the point of the whole 
manoeuvre was that the original Equ. (3.59), when written in terms of hij is 
singular since ct\a — 0. The miracle was that, in the transition from (3.59) to 
(3.61) a factor a is obtained on both sides of (3.61) which can be cancelled 
since a is nonzero outside A by (3.57). 

Thus, taking m sufficiently large and appealing to the theorem of Morrey 
[p5| , we have the 
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Theorem: When M ^ 0, there is a chart in a neighbourhood of A for which 
(a, hij) are analytic. Consequently, from (3.57), Q is also analytic, and so is 

0m = (1-<t)- 3/2 . 

An analogous result can be proved for a suitable set (hij, ft, 4>m, <fis) m 
the stationary case [||, see also [ p9| . The equations one obtains imply in 
particular that the "physical" quantities (hij, <j>Mi <Ps) have an analytic chart 
in a neighbourhood of each point of N and thus entail the "classic" result of 



Muller zum Hagen on the analyticity of stationary vacuum solutions 1 57 . 

By smoothness of (hij , 4>m, ^s) we can define multipole moments for each 
of 4>m,4>s, following (3.46-50). One can show Q that they coincide with 
the quantities E . . . and F . . . in the expansions (3.22-25). (These, in turn, 
coincide with the ones in Thorne j74j], as shown in |24|]). One can now prove 
||, that these moments determine the stationary solution uniquely up to 
isometries. We give a more careful formulation of this result only in the 
static case. 



Theorem: Let there be two static solutions with the same hij\\, the same 
M ^ and the same set of (mass-centered) multipole moments. Then the 
corresponding physical solutions (hij,4>M) are isometric. 

The proof is a not-too-difHcult inductive argument based on (3.61,62), 
(3.68,69) and (3.59,60). 

There remains the question to what degree the multipole moments of 
stationary solutions can be prescribed. It is fairly easy to see, e.g. from 
the asymptotic analysis of Sect. 3.1, that the multipole moments are "al- 
gebraically independent" , i.e. for a given finite number of them, there always 
exists a spacetime having those moments which solves the stationary field 
equations to arbitrary order in 1/r. It is not known what conditions on mo- 
ments for high order have to be imposed in order for the multipole expansion 
to converge. In particular, convergence is not even known when only finitely 
many moments are non-zero. 

There are of course solution-generating techniques to in principle write 
down the general stationary axisymmetric spacetime. To date the only result 
on existence of stationary asymptotically flat solutions without any further 
symmetry is that of Reula |6q | . 

We note, in passing that the above equations lend themselves to an easy 
proof of a result which is often used in black-hole uniqueness theorems (see 
[j3l)). Namely an asymptotically flat, static vacuum solution with M ^ 0, 
which is spatially conformally flat, has to be isometric to the Schwarzschild 
metric near A. To see this, use that now the Cotton tensor of hij is zero. 
Thus, since TZ = 0, the left hand side of Equ. (3.61) vanishes. Contracting 
the r.h. side of (3.61) with (D i a)n jk we find that 

2(jz lj n lj )(Da) 2 = (TZij&u^TVt&a). (3.70) 
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But, by Cauchy-Schwarz, the right-hand side of (3.70) is bounded above by 

which has hence to be zero. Since a can not have critical points near A except 
at A itself, it follows that 

Kii = 0, (3.71) 

whence, from (3.59), D 2 a = 3/2M 2 and thus, in a chart x l for which h%j — Sij 
we have a = \x\ 2 /AM 2 , from which it easily follows that (hij, 4>m) corresponds 
to Schwarzschild with mass M. 



4 Global Rotating Solutions 
4.1 Lindblom's Theorem 

Lindblom showed in his thesis |^| that stationary asymptotically flat dissi- 
pative fluid configuration are axisymmetric. In this section we want to outline 
and discuss this theorem. 

There are three ingrediences of the proof: 

(i) The local fluid field equations imply that the fluid flow is proportional to a 
Killing vector i M provided the divergence of the entropy current vanishes. 

(ii) The Killing field has an extension into the vacuum field of the solution. 

(iii) If the manifold of orbits of the stationary Killing vector £^ is R 3 and 
asymptotically flat, then £ M is linearly independent of i M . The two Killing 
fields commute and there is a linear combination of the two Kiliing fields 
which has fixed points near which it act like a rotation. 



(i) Theorem: Let g^ v , T^ u be a stationary local solutions of the Einstein field 
equations for a one-component fluid with phenomenological heat conduction 
and viscosity laws and vanishing of the divergence of the entropy current. 
Then the fluid flow is proportional to a Killing vector. 

Proof: The energy momentum tensor for a fluid with shear and bulk viscosity 
i s pm (9 and a^ v are the expansion and shear of the fluid; is the heat flow 

M) 

= pu^u v + (p- C9)h^ v - 2ria^ v + q»u v + q v u^ (4.1) 

with 

W u = g »u + u n u v ; q ^ = 0j afij/U n = o _ ( 4 2) 

This implies (a dot denotes the covariant derivative in the direction of the 
fluid flow u M ) 

= -(V„T' U >„ = f> + (p + p)0 - (8 2 - 2^(7^ + V M <f + . (4.3) 



Time-Independent Gravitational Fields 27 



Introducing n, the conserved rest-mass density, and the specific volume v = ^ 
and the specific internal energy u = ^ we can rewrite this, using V /i (n'u A1 ) = 
0, as 

n(u + pv) - (9 2 - 2 W ^a^ + V„g" + g„u" = . (4.4) 

For a one-component fluid we have an equation of state u = u(p,v) and 
consequently there exist scalar functions T(p,v) and s(p, v) with the inter- 
pretation of temperature and specific entropy such that 

du+pdv = Tds . (4.5) 

Hence n(u + pv) = nTs can be used to rewrite (4.4) as 

nTs - (9 2 - 2 W ^a^ + V M g" + g „u" = (4.6) 

or 

ns + T" 1 ^ = T-\Ce 2 + 2na^a^ - g„u") = . (4.7) 
Using again V ^(nu^) = 0we obtain 

V^W+T-V) =T- 1 [(6 2 +2 W ^a^-q^u" +T- 1 WT)} = . (4.8) 

Inserting the phenomenological law of heat conduction 

q^ = -Kh\{T v + Tu v ) (4.9) 

we obtain 

V^nau" + T" V) = T~ 1 ((9 2 + 2 W ^o~^ + kT' 1 ^) = . (4.10) 

The left hand side of this equation is the conserved entropy current V^s^ 
which vanishes according to our assumptions. Hence the positivity of A, ( 
and k implies 9 = o» v = q* = and ii^ = -T _1 T i(U . 

Assume T ^ and consider £ M = T^ 1 ^, the candidate for the Killing 
vector. We have 

V ( „&,) = -T- 2 V iM Tu u) + T~ l W (tl u v) . (4.11) 

The vanishing of 9 = = q» = implies V (ai w„) = -u (m u u) = T^V^T)^), 
hence V( M £„) = 0. 

Now we come to the most complicated part, the extension of the Killing 
vector proportional to the fluid flow from the fluid into the surrounding vac- 
uum region. 

(ii) Conjecture: Let g^ Vl T^ v be a strictly stationary, asymptotically flat 
perfect fluid solution where the matter is a ball of finite extent and the fluid 
flow is proportional to a Killing vector i M . Then f has a unique extension into 
the vacuum region, provided certain differentiability conditions are satisfied 
at the boundary. 
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In Lindblom's original treatment this conjecture was shown to be true 
under the assuption that the outside metric is analytic up to and including 
the boundary S of the fluid. Then one can propagate the Killing vector into 
a neighbourhood of the boundary using the Cauchy Kowalevskaja theorem 
because a Killing vector satisfies a wave type equation. Finally a theorem by 
Nomizu |6l[] can be used to obtain a global Killing vector field. 

One might, however argue, that analyticity up to and including S is too 
strong an assumption. On physical grounds one would like to treat also non- 
analytic equations of state. In this case it is unlikely that the metric is analytic 
in the boundary. 

Finally we show that the new Killing vector i M is actually different from 
the stationary Killing vector £ M . 

(Hi) Theorem: Under the assumption of the above conjecture we have: 

(1) The Killing vectors £ M and < M are linearly independent. 

(2) Both Killing vectors commute. 

(3) There exists a linear combination r/ M = t^ + a^ which has fixed points 
and acts like a rotation with closed orbits. 

Proof: (1) Suppose i M would be linearly dependent of £ M . Then there would 
be a timelike Killing vector, namely , which is asymptotically a translation 
and relative to which the matter does not rotate. Hence, by the Licherowicz 
staticity theorem, spacetime would have to be static. 

(2) As T and m m are invariant objects we have C(T — and C^u u = 
which imply immediately [£,t] = on the support of the matter. To show 
that this is also true outside the matter one can use the analyticity of the 
outer metric up to and including the boundary or one can use a theorem by 
Beig and Chrusciel classifying all possible group action on asymptotically 
flat spacetimes. 

(3) As £ M commutes with i M there is a Killing vector t l on the manifold 
of orbits of £ M . The corresponding group acts in the 2-surface of constant 
pressure, in particular in the boundary, p = 0. As this is topologically S 2 , 
there must be a point where t l vanishes. A Killing vector on a Riemannian 
space with a fixed point acts always as a rotation with closed circular orbits. 
At a point q in spacetime projecting on the fixed point of r, must be 
proportional to £ p and therefore a linear combination if 1 = + with 
constant coefficients vanishing at q exist such that r)^(q) = 0. We have a 
fixed point and because also the timelike direction of £ M is fixed, 77^ acts like 
a rotation and has therefore closed orbits. 

We see that we can obtain the existence of the axis working only on the 
body, provided we know that both Killing vectors are independent. Lind- 
blom |l3| obtains the axis and commutativity of the Killing vectors from 
the asymptotic symmetry group. The key property that the two Killing vec- 
tors are linearly independent is only implied by a global argument and uses 
asymptotic flatness. 
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4.2 Existence of Stationary Rotating Axi-Symmetric Fluid 
Bodies 

Following work by Liapunoff and Poincare, Lichtenstein p2| demonstrated 
at the beginning of this century the existence of rotating fluid bodies in 
Newtonian theory. An account of this almost forgotten work can be found in 
[ [TT) . Using implicit function theorem techniques — as we would say today — 
he shows the existence of solutions near known solutions or approximately 
known solutions: starting with a static fluid ball, he obtains a slowly rotating 
fluid ball; starting from a self gravitating 2-body point particle solution, he 
obtains a solution for two small fluid bodies orbiting their center of mass on 
a circle. Furtheremore, there is a number of exact solutions in Newtonian 
theory: the Maclaurin spheroids, the Jacobi and the Dedekind ellipsoids and 
the Riemann ellipsoids p8[ . 

In Einstein's theory we do not know any stationary exact solution de- 
scribing an extended rotating body. Spacetimes describing such solutions can 
be characterized as follows: Besides a timelike Killing vector £ M there is a 
further symmetry, the axial symmetry generated by rj^, whose orbits are cir- 
cles (Remember that we showed in Sect. 4.1 that such an extra symmetry 
exists on physical grounds) The body is spatially compact and the space- 
time with topology i? 4 is assumed to be asymptotically flat. We assume that 
there is an axis where rf vanishes. Then we can use a result of Carter [jl6| 
which states that under these circumstances the two Killing vectors commute. 
Such spacetimes are called " stationary axisymmetric" , the orbits of the axial 
Killing vector are circles. 

We showed in Sect 2.2 that for stationary axisymmetric perfect fluids 
with an axis and a fluid flow vector contained in the two-surface spanned 
by the two Killing vectors, the two-surface elements orthogonal to the two- 
dimensional group orbit are surface forming ( the group action is orthogonally 
transitive). The same holds in the vacuum region. The property of orthogonal 
transitivity is equivalent to the existence of a discrete isotropy group |70(] . 

To introduce a global coordinate system let us assume that outside the 
2-dimcnsional axis the spacetime is the product of the orbits of the isometry 
group and the orthogonal 2 - surface which we assume to have topology R 2 . 

Using coordinate adapted to the Killing vectors the metric can be written 

as 

ds 2 = g A B(x c )dx A dx B + g w (x c )dt 2 + 2g 0l p(x c )dtd<j) + g rt>rt> (x c )d<j) 2 . (4.12) 

Locally we can always introduce coordinates (x A ) — (r, z) such that g^B 
is conformal to the flat metric in standard coordinates and can therefore write 
the metric as 

ds 2 = e 2k - 2U (dr 2 + dz 2 ) + e~ 2U W 2 dtf - e 2U \dt + Ad<£) 2 . (4.13) 

There is the freedom in (r, z) of an arbitrary conformal transformation 
which is given by the real part of analytic function. 
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The function W 2 is the volume element of the group orbits. As a conse- 
quence of the field equations in vacuum one can locally achieve W = r' such 
that 

ds 2 = e 2k>-2U( dr ,2 + dz ,2j + e -2U r ,2 d(j) 2 _ £ 2U ^ + M ^2 (4 u) 

These coordinates are called Weyl's canonical coordinates. Matters can be 
arranged so that r' = is the axis. Then the coordinates are fixed up to a 
translation in z' . 

It is tempting to try to extend the Weyl coordinates from the outside of 
the body to the interior such that the two-surface orthogonal to the group 
orbit is covered by one (r',z') system with r' = describing the axis and 
W ^ r' in the interior. However, Miiller zum Hagen has demonstrated [ p6[ 
that this is impossible in the case of static spherically symmetric solutions. 
(r' becomes negative inside the body and the axis is reached for p' — > oo.) 
There is no reason to assume that this would be different in the stationary 
case. 

Numerical codes work successfully with a global (r, z) coordinate system 
such that r = is the axis but it is not assumed that one has Weyl's canonical 
coordinate in vacuum. 

For perfect fluids whose velocity is proportional to a constant linear com- 
bination of the two Killing vectors, the case of rigid rotation, V„T Miy = 
becomes particularly simple. (See equation (2.69).) 

= V„T/ = {P + P)\{^r 2 ),» (4.15) 

where / M = / 2 (^ + flrj^) is the four velocity of the fluid. This shows that, 
provided an equation of state p(p) is given, the matter variables p and p can 
be expressed as functions of the quantity / which is determined by the geom- 
etry. This property of rigidly rotating fluids is essential for all the numerical 
schemes as well as for all the attempts to prove existence. 

Various authors have developed codes to calculate numerically stationary, 
axisymmetric rigidly rotating fluid solutions Today this can be done 
with very high presicion by different numerical techniques. These numerical 
solutions are also the basis for investigations of oscillations of rotating stars. 

Schaudt and Pfister try to obtain an existence theorem working in 
the above coordinates adapted to the symmetry. This approach is attractive 
because the field equations become semilinear elliptic. One has, however, to 
control the singularities in the equations on the axis. This is possible and 
two Dirichlct problems have been solved, which give existence of outside, 
asymptotically flat solutions and existence of inner parts of bodies, provided 



appropriate boundary values are given 1 67 . Up to now this was only possible 
for the "outer" and the "inner" problem separately and work is in progress 
which tries to combine the inner and outer solution. 

Let us now turn to the discussion of the only existence theorem for ro- 
tating fluids in Einstein's theory. It is remarkable that the first existence 
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theorem for rotating fluids, proved by Heilig in 1995 |3(J], uses Lichtenstein's 
technique and does not adapt the coordinate to the axial Killing vector to 
avoid difficulties at the axis. 

Let us formulate one particular case of the theorem proved by Heilig . 

Theorem: Let p(p) = Cp 1 be a polytropic equation of state with 1 < 7 < 
6/5. The central density po determines a unique Newtonian static fluid ball 
solution of finite extent. Then there exist a positive constant Qq such that for 
all fl with < fl < Slo a stationary axisymmetric rigidly rotating solution 
with angular velocity fl of the Einstein field equations for a perfect fluid 
exists. The solution is geodesically complete, asymptotically flat with finite 
mass and angular momentum. The matter is of finite extent and has the same 
equation of state and central density as the Newtonian solution. 

The theorem holds also for more general equations of state. It is not clear 
whether the case of positive boundary density may be treated by this method. 

Heilig uses the observation of Jiirgen Ehlers |2l| that it is possible to 
write the field equations as an elliptic system with a parameter A = c -2 - 
interpreted as the velocity of light — such that the equations for , A — > give 
the Newtonian equations and the limit is regular. This can be achieved by a 
particular choice of unknowns for which the field equations are formulated. 

We will describe the structure of Heiligs proof using the equations for- 
mulated in Sect. 2.4 because these are much simpler. We want, however, to 
stress that we expect that Heilig's result could be proved more easily using 
these equations, but this is not certain before all the functional analysis has 
been done properly. 

Let us first adapt the field equations to a rigidly rotating fluid. We write 
the axial Killing vector as in (2.54) 

rj^dx^ = j](cdt + 4>idx l ) + r} i dx i . (4-16) 

The Killing equation in spacetime is equivalent to the equations (2.57-59) on 
the quotient N . For a rigidly rotating perfect fluid with fluid flow vector 
we have 

= + Orin), n = const, u^u" = -c 2 , (4.17) 

where 

f- 2 = e ~^{-e^ + c^flr/) 2 - c- 2 fl 2 i m l . (4.18) 

To obtain the field equation we replace in (2.64-66) ij by c~ x r\ and p by c~ 2 p 
to obtain from (2.46-48) using U — > c~ 2 U 

D 2 U = 47rG[/ 2 (-e^ + c^flrffip + c- 2 p) + 2c~ 2 pe~^} 

+ c - 2 e3f 2 f2 2 {p + c-^rjiTijh^e^ - efwyu" (4.19) 

D%j = 87rGc- 3 e"^/ 2 (-e^ + c^Qi^p + c~ 2 p)f2f]j (4.20) 
Rij = 2c- 4 DiUDjU - 2e%QikQj k + h l3 e^uj k iuj a 
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+ 87rG<r 4 h2pe~^ y - + f 2 2 (p + c _2 p)»R% 

- hijf 2 n 2 (p + c- 2 p) mr]m h lm ] (4.21) 

The above field equations have to be supplemented by the Killing equations 
(2.57-2.59). The equations of motion are 

V^T"" = <=> (c 2 p + p)f~ 1 D i f = D iP (4.22) 

For c-toowe have from (4.18) that f 2 = l which implies by (4.20) that 
D l u>ij = 0. The staticity theorem now gives that oXy = 0. Using (2.58) this 

_1U 

implies Di{e "^rj) = 0. The vanishing of ?y on the axis implies r\ = 0. Using 
all this the field equations reduce to 

Rrj = , D 2 U = AnGp (4.23) 

Therefore the metric on N is flat. Using Zim c _>oo(c 2 -Di/) = —Di(U~^Q 2 r]iri l ) 
the equation of motion become the Newtonian equation 

-pA(f/-^W)=Ap (4.24) 

equation (rjiT] 1 = x 2 + y 2 in Cartesian coordinates). 

As discussed in Sect. 2.5 for the static case, the field equations become 
again a quasilinear elliptic system for U, Z' L ° , tpt in harmonic coordinates 
(V^Vi = <^=> D i (p i = 0, V^Vx* = D 2 x l = 0). Namely, the condi- 
tion that the time function is harmonic turns the left-hand side of Equ.(4.20) 
into an elliptic operator acting on tpi. Harmonicity of x % has the same effect 



on the left hand side of Equ.(4.21). Theorem 4.1 of Heilig 30 can be adapted 
to show that for small A, £2 a solution of the reduced field equation exists near 
the Newtonian solution. Such a solution satisfies only the harmonicity condi- 
tion if the equation of motion holds. So, this has to be solved simultaneously. 
This is possible because given a equation of state (4.22) can be integrated 
such that the matter quantities can be expressed in terms of the geometrical 
quantity /. Therefore the following iteration procedure is well defined: begin 
with a Newtonian solution U°,p°; choose some A, f2 and use p°,p°, U° as a 
source in the field equations in harmonic coordinates to obtain U , Z 3 , ip X i. 
Calculate / from U 1 ,Z U - i,ip l i,\,Q and determine p 1 from the equation of 
motion. Then one solves again the field equation with the new source and 
so on. Heilig has shown that for sufficiently small A and fl such an iteration 
converges in his variables. It should also converge in the variables used here. 

It is remarkable that we have used rjidx 1 = xdy — ydx as a given field. At 
the end one has to check that the solution is axisymmetric and satisfies the 
harmonicity condition. 

Note that only for A = cr 2 with some fixed value of the velocity of light 
in some units the above field equations are Einstein's equations. It is however 
possible to reinterpret solutions with any A as solutions of Einstein's equation 
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expressed in different units |3Cj • With this interpretation the theorem above 
demonstates the existence of slowly — the theorem does not control the range 
of u — rotating fluid configurations. 



4.3 The Neugebauer Meinel Disk 

The only known global solution describing a rotating object in Einstein's 
theory, is the relativistic analogue of the rigidly rotating Maclaurin disk in 
Newton's theory pj. 

An axisymmetric surface density distribution (in cylindrical coordinates 

o-(r) = ctqa 1 - T — , 0<r<a, (4.25) 
V r o 

generates a gravitational potential <P(r, z), which is determined by the Poisson 
integral from a. At the disk the potential is 

*(r, 0) = -Q 2 r 2 + const , < r < a , Q 2 = . (4.26) 

2 r^ 

Outside the disk the potential can be expressed, for example, in terms of 
integrals over Bessel functions. 

The centifugal force acting on rigidly rotating particles balances the grav- 
itational force in the disk. Therefore, we can interpret the density distribu- 
tion as formed by self gravitating, rigidly rotating dust. The two parame- 
ters (To and ro determine such disks uniquely. The total mass of the disk is 
M = ^iraor 2 . 

Neugebauer and Meinel found the relativistic analog of these disks |5j| . 

There is a well known formalism available in General Relativity to de- 
scribe matter surface distributions |nj . In the particular case of a reflection 
symmetric disk, we have to find solutions of the stationary vacuum field 
equations, defined outside the disk such that the difference of the normal 
derivatives of the metric at of the disk have a certain algebraic structure ]3l| . 

The general stationary axisymmetric metric can be parametrized as 

ds 2 = e- 2U [e 2k (dr 2 + dz 2 ) + r 2 dtf] - e 2U (dt + ad<j>) 2 . (4.27) 

The metric coefficients U, k and a depend only on r,z; the vector fields 
^d/dx^ = d/dt and rfdjdx^ = d/d<j> are Killing vector fields. We assume 
that the orbits of the axial Killing vector are circles; r — is the axis. 

Let us assume that the disk is located at z = 0, < r < ro. A rigidly 
rotating flow forming the disk is described by a vector field (which is defined 
at the disk) 

= e~ v {e + Qr}") , u^u^ = -1 , (4.28) 

where fl is constant. The definition of a dust disk implies that the metric 
is continuous across the disk and that t^ v = au^u w , where a is the surface 
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density, satisfies t^" - v — with respect to the Levi Civita connection of the 
metric induced on the disk. As = £ M + Q<nf is a Killing vector, it holds 
iA p = 0, <7 4l v^ — 0, V^v^ = and we obtain 

t" u . v = (ae- 2V , v = ae- 2V v tt . v v v . (4.29) 

Finally e 2V = g iiV v ll v v implies e 2V 2V 7 = 2g liv v fl v 1 ' ;7 = 2v v v 1 - u and we see 
that V must be constant on a disk formed of rigidly rotating dust, V — Vq. 
It is natural to introduce comoving coordinates 

t' = t; 4>' = <t>- Qt , f' + Mr/'* , if' = if , ^ = e~ V St'. (4.30) 

The vacuum field equations can be expressed in terms of the following quan- 
tities: 



< - = = e 2V , a' = - e - 2U 'n^' , U'(r, ( f ) ,z = 0) = V Q = const 



(4.31) 

and b'(r,z) determined by 



a\ r = re- 4U 'b[ z , a[ z = ~rer w ' b\ r . (4.32) 

Using the Ernst potential /' = e 2U + ib' the key field equation is the semi- 
linear elliptic Ernst equation J37J 

Mf')(f',rr + /,« + \f,r) = C + /* ■ (4-33) 

For a solution of the Ernst equation the integrability condition of (4.32) is 
satisfied and one can solve for a' . The remaining metric coefficient k' follows 
from the equations 



u' 2 -U' 2 + \e^\b 2 -b 2 z ) 



k' = 2r 



U[ r U[ z + \e^ u '{b[ r b[ z ) 



(4.34) 

whose integrability condition is again satisfied for solutions of the Ernst equa- 
tions. 

We can perform an integral in the (r — z) -plane around the disk of the 
integrability condition of (4.32), namely 

(r-VVrXr + (r-V^Vj,* = , (4.35) 

which can be replaced by a surface integral. As we assume that the tangential 
derivatives of the metric are continuous at the disk, we obtain at the disk 

< z Uo+ = <J*=o- • (4.36) 
On the other hand reflection symmetry at the disk implies 

<rUo+ = ~<z U=o- (4-37) 
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on the disk, hence, 

<,U=o+ = <U=o- = , (4.38) 

which by (4.32) implies b' = const on the disk. 

Now it is easy to calculate the second fundamental form k cd = g c d,z 
of the disk z — ( c,d,... = (t,r, </>)), because we have at the disk that 
a' z = k' z = 0, as a consequence of (4. 38), (4. 35) and (4.34). We find 

k rr = -2U' z g rr (4.39) 

= -2U' <z (a' 2 e 2U ' + e- 2C V) (4.40) 

h't' = W[ z g vtl (4.41) 

k V p=2U[ z gw (4.42) 

Now we can check the condition for a disk of dust , namely 

+ k cd - ~k cd = 2 + k cd = -87r(r cd - -Qcdji) = -%n(<ju c u d + ^ag cd ) , (4.43) 

which, in the primed coordinates (because of = 5?, ) , reads 

k c 'd' = -8Ka(g c > t ig d it> + g c 'd') ■ (AAA) 

Because of the form of the metric (4.27) in primed coordinates and by (4.39- 
42), this is satisfied if we define the surface density by 

a = ±-U[ z . (4.45) 

Z7T 

Thus we have shown that a rigidly rotating disk of dust is determined by 
a solution of the Ernst equation which satisfied at the disk U' = const and 
b' = const. Outside the disk the solutions of the Ernst equation must be 
regular. For a well-posed elliptic problem we need furthermore asymptotic 
flatness at infinity and regularity conditions at the axis. 

In Newton's theory there is a 2-parameter family of disks (4.25), (4.26). 
If we use the 2-parameter group of similarity transformations — or dimen- 
sionless quantities — we can assume ro = 1 and ao — I and we have just one 
disk. 

Because of the appearence of the velocity of light there is only a 1- 
parameter group of similarity transformations in Einstein's theory. Hence, 
after we put rg = 1, we expect a 1-parameter family of disk solutions. 

The investigations of Neugebauer and Meinel suggest that 

H = 2uj 2 rle- 2Va (4.46) 

is an appropriate parameter. 

Neugebauer and Meinel prove by the so called inverse scattering method 
of soliton physics (compare the contribution of Maison in this volume) that 
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the boundary value problem for /' has a unique global solution provided Vq 
and lo are such that /i < [i cr it = 4.629 . . .. 

The solution /' can be expressed in terms of hyperelliptic theta functions 
[ |60| . The remaining metric coefficients al and k' are determined by integration 
from the equations (4.32) and (4.34). 

If we put the velocity of light, c, in the appropriate places we obtain the 
MacLaurin disk as a Newtonian limit. 



Further properties of these disks are discussed in 




Many global stationary solutions with disk sources may be constructed 
from known stationary vacuum solutions by "cutting out" a region containing 
singularities an making appropriate identifications. This method was first 
used by Bicak and Lcdvinka || to produce physically plausible sources for 
the Kerr metric with arbitrary values of the parameters a,M. These disks 
are made of two streams of particles circulating in opposite directions with 
differential velocities. They are extending to infinity but have finite mass. 
See Sect. 6 of the article by Bicak in this volume, where this procedure is 
related back to the "method of images" in Newtonian galactic dynamics. In 
the static case these methods yield an infinite number of such static disk 
solutions. Solutions corresponding to stationary counterrotating dust disks 
of finite extent have been constructed by Klein and Richter |35| . 

5 Global Nonrotating Solutions 
5.1 Elastic Static Bodies 

No doubt, Einstein's theory should allow for the description of static, solid 
bodies. It is useful to make the following distinction: 

(i) small bodies, whose shape is not dominated by gravitational forces, 
like a piece of sugar or an iron ball. If we ignore gravity, the structure of 
the body is determined by the laws of quantum mechanics. This is true in a 
Galilei invariant formulation as well as in a special rclativistic one. Linear and 
nonlinear elasticity theory describes the deformation of such a configuration 
under external forces. 

Suppose we now want to add the gravitational field. This is straightfor- 
ward for linear elasticity in Newtonian theory; we just have to insert the 
gravitational field calculated from the Poisson integral as an external force 
into the equations of elasticity. 

To pass from special relativity to Einstein's theory is more complicated. 
Now the deformed configuration has to satisfy Einstein's field equations, and 
the elasticity equations are a consequence of the latter! 

(ii) bodies like stars whose shape is dominated by gravitational forces. 
There a relaxed state does not really exist and one has to modify the de- 
sciption of elasticity. This holds in Newtonian theory as well as in Einstein' 
theory. 
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Elastostatics can be described in Einstein's theory as follows |i7|| . The col- 
lection of particles which form the body is described by the three-dimensional 
"body manifold" B. The essential dynamical variable is a map <P : M A — » B. 
such that <P^ 1 (y l ) is the world line of the particle in spacetime labeled by 
y l in B. In the static case the world lines of the particles are the integral 
curves of the Killing vector, and we can consider <P as a 1-1 map N — ► B. 
We assume that we have given on B a Riemannian metric Rij. Its physical 
interpretation may be different: for small bodies it describes a relaxed state; 
for big bodies which go never into a relaxed state, it could be an "isotropic 
state of minimal energy" . 

We need now information about the energy momentum tensor of the 
material. Let 

= pu»u v + p» v , p» v u u = (5.1) 

be such that the stress tensor p^ v has only spatial components and can be 
considered as a tensor on TV. We can now define 

e ij : = 2^ hij ~ ®* Ri i) ( 5 - 2 ) 

the "Lagrangian strain tensor" . In the Hookian approximation of elastcity 
one assumes that one has given on the body B a tensor field K v such that 
after moving this object with <1> into the space N' one can define 

P = Po + ^K ijki e i:j eM (5.3) 

p V = -K^ M e kt (5.4) 

as the energy and stresses of the body B in 3-space with the metric hij . 

With this energy momentum tensor we consider Einstein's field equations 
as differential equations for the spacetime metric and the map No general 
existence theorem is available for this problem. The only case treated so far 
is the spherically symmetic one Jfj^] . 

To get some feeling for these equation let us consider some further ide- 
alisation. For small deformations we can linearize := h{hij — <?*K.y) as 
follows: Suppose that \P € is a 1-parameter family of diffeomorphism N 1 — > N 1 
such that e = is the identity and is some diffeomorphism N' — > B. Now 
we assume that ?P r e # _1 o defines our deformed body and calculate the stress 
tensor to first order in e If we define = <?o*Sy we obtain 

e* = \{hij - k% + £ x <.) = \{h l3 - 4 + D\ a k K% h ) (5.5) 

Here the vector field \ l is defined by the linearization of ^ on N' . We see 
that p l J -j — leads to second order differential equations for % a . 

Consider first the case of special relativity which coincides with Galilei 
invariant classical mechanics in the static situation. Then we have h a p — rj a p 
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and hij = 6ij. We choose to be the identity may which discribes the 
relaxed body in spacetime. With R{j = Sij we obtain 

This gives implies the equations of classical, linearized elastostatics . 

/'''., A''" M.jj (5.7) 

With the appropriate symmetry and positivity conditions on the equa- 

tions are elliptic and solutions exist for various boundary conditions. 

Next we want to calculate the deformation of a small elastic body by 
its own gravitational field. The relaxed state is determined by solid state 
physics as above. To switch on gravity we assume that we have families 
g„u = r)„u + Gg 1 ^ + G 2 gl„ . . . and = T° v + GT^ + G 2 T^ . . . satisfying 
the field equations. 

At order G° we obtain the trivial solution if there are no forces at the 
body , the density pg 1S constant and the stresses vanish, i.e. \a — 0- The field 
equation in order G are obtained from the equations in section 2.4 with an 
energy momentum tensor T® v which has only a term p® because the stresses 
vanish. We obtain U 1 as a solution of the Poisson equation with the source 
source p$. The metric hjj remains fiat in this order. The expansion of the 
equation of motion in G gives to first order 

P lij ,i = -K ijH xl tj = -elU\ a , AXJ 1 = AnGp° Q (5.8) 

Hence we obtain classical elastostatics with the force deforming the body 
being the gravitational force. 

One might try to obtain an existence theorem for small self gravitating 
elastic bodies in Einstein's theory by an implicit function theorem argument 
similarly as in the case of a rigidly rotating body (section 4.2). 

5.2 Are Perfect Fluids 0(3)-Symmetric? 

It is intuitively "obvious" that a static, in particular non-rotating, ball of 
perfect fluid, due to the absence of shear stresses should have spherical sym- 
metry, and in particular the gravitational field in its exterior should be the 
one described by the Schwarzschild spacetime. This result, in its most general 
form, is still open in G.R. (The Newtonian case was settled in Lichtenstein fl42| 
and Carleman jl5), see also Lindblom Q|.) Rather, one has today a theorem 
which is essentially a uniqueness result in the spirit of black hole uniqueness 
theorems. An earlier result due to Kiinzle and Savage Q states that, near a 
spherical solution, there is no asphcrical one with the same equation of state 
and the same mass. 

The following result, due to Beig and Simon B , is a refinement of previous 
work by Masood-ul-Alam |52| , see also the review of Lindblom pjj . 
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Theorem : Let us have a static, asymptotically flat, spherically symmetric 
solution to the Einstein equations with a perfect fluid and barotropic equation 
of state p = p(p). (This solution is called reference spherical solution.) Let 
there be given another static, asymptotically flat solution with the same 
equation of state and the same value V\qs of the Killing vector norm on 
the surface dS of the star. Let further p{p) satisfy the differential inequality 
/ < 0, specified later. Then these two spacetimes are isometric, in particular 
the second one is also 0(3)-symmetric. 

The condition stipulating the existence of a spherical reference solution 
was disposed of by Lindblom and Masood-ul-Alam . The condition on the 
matter, besides the one stating that p > 0, p > and dp/ dp > 0, is that 

1 dn 

I := -k 2 + 2k+( P + p)— <0, (5.9) 
5 dp 

where k := f+§^ One can check that it is for example satisfied for the 
equation of state of a relativistic ideal Fermi gas at zero temperature, but 
only up to densities of order 10 15 gcm -3 , which is roughly the critical density 
where gravitational instability sets in. It is known from numerical results [^9| 
that beyond that limit the uniqueness statement of the above theorem will 
fail. One believes however, that sphericity will still hold. 

We will here confine ourselves to an outline of the proof to the case of the 
special equation of state given by Q 

p(p)^lp e/5 (p 1 o /5 ~P 1/5 r 1 (P0 = const >0) (5.10) 

which is a relativistic generalization of the equation for a polytrope of index 
5 in the Newtonian theory . The expression in (5.10) satisfies 1 = 0. The 
reference spherical solution in this case is known explicitly [Q. It has the 
curious property that it is asymptotically flat, but the fluid extends to spatial 
infinity. 

Introducing the variable v = (— V) 1 ^ 2 , the static field equations for a 
perfect fluid with energy momentum tensor 

Tpty = {p + p)u^u v + pg^ (5-11) 

with w M = v^ 1 ^ read 

D 2 v = 4irv{p + 3p) (5.12) 

Tlij = v~ x DiDjV + 4vr(p - p)hij (5.13) 

The asymptotic conditions (2.102,104) imply that v — ► 1 at infinity. ^From 
the maximum principle for elliptic equations it follows that < v < 1 in N. 
Since the surface of the star is at infinity for the Buchdahl solution, the v\qs, 
which is always equal to one in that case, has to be replaced by the total 
mass M (see Sect. 4.1). 
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Applying the contracted Bianchi identity to (5.12,13), there follows 



Dip = -v 1 (p + p)D i v. 
Thus p and p are both functions of v and 

-v- x {p+p). 
Define the Cotton tensor of hij B^k by 



dp 
dv 



[k 



With the definition 



W = (D lV )(D l v) 



(5.14) 
(5.15) 

(5.16) 
(5.17) 



the equations (5.12,13) now imply (see Lindblom |46 ) that 



D 2 W = -v i W- 1 B ijk B ijk + v- 1 {D l v){D l W) + %nv{D l v){D t p) + -f^ffffJtAR') 



- 8nW(p + p) + 16n 2 v 2 {p + 3p) 2 - Anv{p + 3p)W- 1 (D l v)(D l W). 

In the spherically symmetric case W = Wo has to be of the form Wo — Wq(v). 
The ODE resulting in that case from (5.18), has, for the equation of state 
(5.10), an explicit solution given by 



(5.18) 



Wo 



,^4 



16M 2 



3 



(5.19) 



We assume that a = poM 2 > 1. The function Wo is defined for v G [0, 1]. 
It is positive for v £ (v c , 1), with v c = (yfa — l)/( v / a + 1) and Wo(v c ) = 0, 
^(l) = 0. Thus Wo satisfies the correct boundary condition at the central 
value v c of v and at infinity. 

We now define, for the given solution (v, hij), the scalar function 



W -Wo 



and the conformally rescaled metric 



hi 



(W - Wo) 



(5.20) 



(5.21) 



(The constant M occurring in Wo is taken to be the mass of the given so- 
lution.) In the asymptotically flat, vacuum case discussed in Sect. 3.2 one 
finds that the metric hij extends smoothly to the manifold N — N U {A}, 
with A the point at infinity. This is also true for the Buchdahl solution, and 
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we assume it to be true for the given, a priori non-spherical one. After some 
calculations we find that 

D 2 (W - W ) = -^B ljk B %3k + ^W^&iW - W )Di(W - W ). (5.22) 

Since W, W also extend smoothly to N, the function W — W satisfies the 
elliptic equation with nonnegative right-hand side on the compact manifold 
N. After integrating (5.22) over N (or by the maximum principle) it follows 
that Bijk is zero and 

W = W (v). (5.23) 

It then follows from that the given model is isometric to the Buchdahl 
solution with the same po and the same M. 



5.3 Spherically Symmetric, Static Perfect Fluid Solutions 

The metric for a static spherically symmetric spacetime can be written 

ds 2 = ~c 2 e v{r) dt 2 + e x{r) dr 2 + r 2 (d9 2 + sin 2 9 d<t> 2 ) . (5.24) 

For a derivation see (2^j7^|. Here c is a constant which plays the role of the 
speed of light. In Appendix B of |^8| it is also demonstrated that the r 2 in 
front of the sphere metric is no loss of generality for a static perfect fluid 
with positive mass density and pressure. Hence it is impossible to have two 
centers or two infinities. The field equations for a perfect fluid are 

8ttGc- V 2 = e~ x (rX' - 1) + 1 (5.25) 
87rGc~V 2 = e- x {rv' + 1) - 1 (5.26) 

8ttGc~ 4 P = X -e- x (v" + X -v 12 + r" V - A') - ^v'X^j (5.27) 

A prime denotes a derivative with respect to r. We have written —c 2 p for the 
timelike eigenvector of the energy-momentum tensor to make the compari- 
son with the Newtonian equations easier. The field equation imply 'energy- 
momentum conservation', which is a single equation for a static perfect fluid 

2p' = -v'{p + c 2 p) . (5.28) 

The first exact solution of these equation was alredy found in 1918 by Karl 
Schwarzschild, the solution with constant density |37|. We have three in- 
dependent ordinary differential equations for for four functions. Hence, one 
function can be specified freely. The most physical case is to prescribe an 
equation of state p = p(p). Equation (5.25) can easily be integrated: 

er x = 1 - ^r-i f r 2 p(r)dr + const . (5.29) 
& r J 
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As we only are interested in solutions with a regular center of spherical sym- 
metry we define A as follows 

e -A = 1 _ f r 2 p(r)dr . (5.30) 

c z r J 

The usual definition of the 'mass up to r\ namely 

m(r) = 4tt / r 2 p{s)ds (5.31) 



gives 

e - = l-^l. (5.32) 

It is also useful to introduce the following quantity which is related to the 
'mean density up to r' 

w(r) = r~ 3 m(r) . (5.33) 

Then (5.32) becomes 

e" A = 1 - ^-r 2 w . (5.34) 

Various forms of the equations (5.25-28) will be used. Equations (5.25), (5.26) 
and (5.28) contain all the information. If we eliminate v' then (5.26) and 
(5.28) imply the Tolman - Oppenheimer - Volkoff equation fff§ 

( 2G 2 \ _1 (^p , \ (P 
p =-Gr[ l-—rw) -r 



w Mj + p (5.35) 



If an equation of state is given we can integrate (5.28) 

/■pO) 2dp 

v(r) = „ , , + constant (5.36) 

J P0 P + c 2 p{p) 

In this formula po denotes the central pressure. If we add the definition of 
w then (5.35) and (5.33) form an integro-differential system. Differentiating 
(5.33) we obtain 

w' = -Uirp-3w) (5.37) 
r 

In |65| the following theorem is proved. 



Theorem: Let an equation of state p(p) be given such that p is defined for 
p > 0, non-negative and continuous for p > 0, C°° for p > and suppose 
that dp I dp > for p > 0. 

Then there exists for any value of the central density po a unique inexten- 
sible, static, spherically symmetric solution of Einstein's field equation with 
a perfect fluid source and equation of state p(p). The matter either has finite 
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extent, in which case a unique Schwarzschild solution is joined on as an ex- 
terior field, or the matter occupies the whole space, with p tending to as r 
tends to infinity. 

There are two parts of the proof. The equations (5.35) and (5.37) form a 
system of ordinary differential equations for p(r),w(r). However, the system 
is singular at r = and the first step is to demonstrate that for each value 
of the central density there is a unique solution such that the spacetime is 
regular at the center. This is shown in j65| or in (jso). This solution defines a 
neighborhood of a regular center and can be extended as long as (1 — ^-r 2 w) 
remains positive. This can be seen as follows. 

Introduce the variables first used by Buchdahl [[l3) 

y 2 = 1 - ^r 2 w , C = <? /2 , x = r 2 (5.38) 

Rewriting the equations in these variables and eliminating p in 5.26) and 
(5.28) gives an equation which is linear in £ and w, 

(1 -^-xw)Cxx - %C,,x{w + xw, x ),x - 7T^W, X (, = (5.39) 

or 

0/C,;., tj~t ( 5 - 4 °) 

Let < x < xo be an intervall such that y 2 — (1 — ^sro > and p > 0. As 
the density does not increase outwards we have w jX < 0. Therefore 

(»C*),x < (5.41) 
The equation (5.26) can be rewritten as 

y ^ = y^ W + $ p) ■ (5 ' 42) 
^From (5.41) and (5.42) we obtain the inequality 

w + Airp/c 2 

V > — — : ■ (5.43) 

w + 4:np Q /c z 

Hence, we see that y cannot vanish before p. 

Suppose p{xb) = 0. Then we call the corresponding the radius of the 
star. The Schwarzschild solution is given in the form e~ A = e v = 1 — A/r 
for some constant A. Hence, we determine a unique exterior field by the 
condition A = (^r)m(r-fc). In this way the matter solution and the outside 
solution are joined only in a C°-fashion because the boundary density may 
be non-zero. If we introduce Gauss coordinates relative to the hypersurface 
p = the metric is C 1 . It is obvious that this metric cannot be extended 
because the area of the group orbits r = constant grows from to infinity. 
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Let us now consider the second possibility that p(x) > for all x. Because 
p{x) is monotonically decreasing for x — > oo, lim^^oo p(x) — p^ exists. This 
implies that p' tends to for x — > oo. Since y < 1, (5.35) then implies p^ = 
and hence, using the equation of state, that p — ► as x — > oo. As before the 
spacetime is not extensible. 

This completes our outline of the proof. It shows in particular that for for 
p(p) with p(0) = pi, > the radius of the star has to be finite. 

There are various exact global solutions known. (For a useful list of such 
solutions including a discussion of their physical acceptability has been given 
by Delgaty and Lake @].) For the 1-parameter family of equations of state 
given by Equ. (5.10) the whole 1-parameter family of solutions is known. A 
2-parameter family of equations of state of interest for the issue of section 
3.1 is investigated by Simon |fr2f ; all the corresponding exact solutions are 
given. 

There are some conditions on the equation of state known, which allow 
to decide whether the radius of the star is finite or infinite in the case of 
vanishing boundary density pb- In |65| it is shown that the radius of the star 
is finite if dp/ p(p) 2 is finite. Conversely, JJ° dp/(p(p)c~ 2 p < oo implies 
that the matter distribution is infinitely extended. Both conditions depend 
only on the behaviour of the equation of state near the boundary p = 0. 
Makino |30[ gives conditions for a finite radius in cases which are not covered 
by the above. He shows in particular, that for polytropic equations of state, 
p = const. p 1 with 4/3 < 7 < 2 the radius is finite. 

For finite distributions " Buchdahl's inequality" holds (l3| . 

Theorem: For finite distributions with non-negative density and a monotonic 
equations of state there holds 



!-——>-, (5.44) 
c r& 9 



Proof: To obtain the inequality one compares the solution with a solution of 
constant density p, an interior Schwarzschild solution. Equ. (5.40) implies for 
this solution (written with an overbar) that 

(yC,x),x = ==>• y( :X = a = constant (5.45) 

We normalize C by the condition that at the boundary we have Cb — yb- Then 
we find a if we rewrite (5.26) in the new varables 

8nG a 2C.X 2G 

— 2~P = 4y — T w (5.46) 

and evaluate it at the boundary as a = ^§rW. 
Then (5.45) can be integrated with the result 



C(*) = Ml + 2C(0) + yi-|^}) (5.47) 
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Now (5.41) implies 

yC.x > (y(,x)b = vC,x (5.48) 

As y > y we obtain 



C*>Cx = Hl + 2C(0) + yi-^!B«;J (5.49) 



As C is positive we obtain at the boundary 



Vb>-\w>+\ (5-50) 



which is (5.44). 

Buchdahl's inequality show that one can pack only a certain mass into a 
given fixed radius. The physical reason is that the pressure is also a source 
of the gravitational field. In Newton's theory there are constant density balls 
with a fixed radius for arbitrary density. In Einstein's theory the central 
pressure diverges if the density approaches some maximum value. 

In pi] an analogue of Buchdahls inequality is derived for distributions in 
which the the density is only assumed to be positive. There holds 1— > 
0. 

Another important topic are bounds on the total mass of the system. 
Suppose we know the equation of state only for p < p . Then we can estimate 
the mass and radius of a core in which the density is greater po as follows: 
Clearly, m(r ) > ^-po(r ) 3 ; because of y > we have also m(r ) < §qTq. 
Hence the possible cores occupy a compact part of the m(ro)-ro _ plane. 
Taking intial values from this part one can numerically integrate outwards 
using the known equation of state, until the pressure vanishes. This was done 



in 1 27 for po = 5.1 x 10 1 g/cm 3 and with a certain realistic equation of 
state for smaller densities. All configurations had a total mass smaller then 
5M . It is quite remarkable the the knowledge of the equation of state for a 
finite density range allows to show such a bound on the total mass, assuming 
nothing but the monotonicity of the equation of state in the unknown density 
range. This is not possible in Newton's theory. 

In the special case of bodies with a sharp edge, i.e p^ > 0, we can combine 
the Buchdahl inequality (5.44) with the estimate M = m(rb) > Anpbr 3 to 
obtain the mass bound 

-(!)W 2 - 

Let us finally compare with Newton's theory. In (5.35) it is almost obvious 
that for c — > oo one obtains the Newtonian equation for the pressure. The 
relativistic corrections show how "the pressure enters in the active and passive 
gravitational mass" . The first factor describes an effect of the geometry. Static 
fluid ball are the simplest examples of families of relativistic solutions with a 



Newtonian limit |21 
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5.4 Spherically Symmetric, Static Einstein Vlasov Solutions 

In recent years existence and further properties of solutions of Einstein's 
field equations for a collisionless gas have been shown Q. The Vlasov- 
Einstein- system determines the spacetime metric and the distribution func- 
tion /(a^jj/ 1 ) describing the particles. 



T"" = / P ll p v \g\ 1/2 ^ (5.52) 
= 8irT^. 

In the static spherically symmetric case and for the metric (5.24), these equa- 
tions reduce to (r= are the spatial frame components of p a ) 

—L=d x ,f -VT+ «V-cW = (5.53) 
VI + v 2 r 

8TrGc- 2 pr 2 = e- x {r\' - 1) + 1 (5.54) 
8ttGctV 2 = e~ x (rv' + 1) - 1 (5.55) 



where 



p(x) = p(r) = / f{x\v i )^l + v 2 dv , (5.56) 

J R 3 

p(x)=p(r) = / fixW) (—) -jP=^ ■ (5-57) 

The distribution function is assumed to be spherically symmetric. 

Rein and Rendall |i4j show the existence of asymptotically flat solutions, 
regular at the center, with finite total mass and finite extension of the mat- 
ter and isotropic pressure. It is also possible to construct solutions with 
anisotropic pressure; Furthermore shells of finite extent of matter around 
a regular center or a black hole can be constructed (63j . 
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